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Neutrino emissivities due to direct Urea processes of several spin-one color-superconducting phases 
of dense quark matter are calculated. In particular, the role of anisotropics and nodes of the gap 
functions is analyzed. Results for the specific heat as well as for the cooling rates of the color-spin- 
locked, planar, polar, and A phases are presented and consequences for the physics of neutron stars 
are briefly discussed. Furthermore, it is shown that the A phase exhibits a helicity order, giving rise 
■ to a reflection asymmetry in the neutrino emissivity. 
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Exotic states of matter are expected to exist in the central regions of compact stars. The baryon density in these 
systems is likely to exceed several times the nuclear saturation density, po ~ 0.15 fm -3 . The exact nature of such dense 
matter, however, is yet unknown. It was suggested already 30 years ago Q that it might be deconfined quark matter. 
\ Since the temperatures in neutron stars are sufficiently low, this matter is likely to be in a color-superconducting 
state. (For reviews on color superconductivity see Refs. 0.) There is little doubt that dense baryonic matter at 
asymptotically high quark chemical potential fi is a color superconductor In this case, the ground state is the 
color-flavor locked (CFL) phase (for studies of the CFL phase in QCD at asymptotic density, see also Refs. 0,0)- 
At densities existing in stars, however, this phase may not necessarily be realized. The main reason for the potential 
breakdown of the CFL phase is a relatively large difference between the masses of the strange quark, m s , and the 
masses of the up and down quarks, m u ~ mj, which is negligible only at large densities. This difference, together 
with the requirements of (3 equilibrium and electric and color charge neutrality, gives rise to a mismatch between 
the Fermi momenta of the quarks that form Cooper pairs Q. Therefore, the conventional BCS pairing which is 
the underlying mechanism of the CFL state, becomes questionable and the search for the true ground state of quark 
matter in compact stars has to include different, unconventional superconducting states |9l Hoi fllL Il2| . 

Determining the ground state of QCD at moderate densities from first principles, i.e., within the framework of 
QCD, is not possible at present. It is natural therefore to make use of astrophysical observations to test the presence 
of various suggested phases. The goal of this paper is to investigate color superconductors in which quarks of the 
same flavor form Cooper pairs and to study their effect on specific observablcs. 

Single-flavor Cooper pairing is the simplest possibility for neutral quark matter. Contrary to other unconventional 
pairing mechanisms, such as the gapless 2SC phase 0, the gapless CFL phase the CFL phase with additional 
meson condensates |llj . or the crystalline phases [l2j |. it is allowed for arbitrarily large mismatches between the 
Fermi momenta of different quark flavors. Single- flavor pairing in the color anti-triplet channel is possible only in 
the symmetric spin-one channel 0, El 0, 0, Gj| El ■ This is the consequence of the Pauli principle which requires 
the wave function of the Cooper pair to be antisymmetric under the exchange of the constituent quarks. This is in 
contrast to pairing of different flavors where the antisymmetric spin-zero channel is allowed. 

Furnishing triplet representations with respect to both color and spin groups, the order parameter in a spin-one 
color superconductor is given by a complex 3x3 matrix. This is similar to supcrfluid 3 He, where condensation occurs 
in spin and angular momentum triplets |l8j| . In both systems, the matrix structure of the order parameter gives rise to 
several possible phases. In 3 Hc the observed phases are the so-called A, B, and A\ phases. In this paper, we consider 
the following four main spin-one color-superconducting phases of quark matter: the color-spin locked (CSL), planar, 
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polar, and A phases, proposed in Refs. [Til HM IT^. ITtI ITgj. The A and CSL phases are analogues of the A and B 
phases in 3 He, respectively. 

In contrast to the nonrclativistic case of 3 He, spin-one color-superconductors may involve pairing of quarks of the 
same as well as of opposite chiralities. In this paper, we focus on the so-called "transverse" phases, in which exclusively 
quarks of opposite chiralities pair. Theoretical studies at asymptotically large densities suggest that these phases are 
preferred (lllT7j . 

As in the case of 3 He, the gap functions of most of the quark phases considered here are anisotropic in momentum 
space. Only the CSL phase is isotropic. The gap in the polar phase vanishes at the south and north poles of the 
Fermi sphere, whereas the gap in the planar phase is anisotropic but nonzero in any direction of the quasiparticlc 
momentum. The A phase is special in the sense that it has two gapped quasiparticlc modes with different angular 
structures, one of which has two point nodes. 

One should expect that the anisotropies and especially the nodes affect the physical properties of the corresponding 
quark phases. The low energy excitations around the nodes give important contributions to various thermodynamical 
and transport properties, e.g., the specific heat, the neutrino emissivity, the viscosity, the heat and electrical con- 
ductivity, etc. In application to compact stars, the specific heat and the neutrino emissivity determine the cooling 
behavior during the first 10 5 - 10 6 years of the stellar evolution. In this paper, we compute these two quantities for 
the four mentioned spin-one color-superconducting phases and deduce the resulting cooling rates. In addition, we will 
discuss the reason why the distribution of neutrino emission from the A phase breaks reflection symmetry in position 
space. 

The paper is organized as follows. The general formalism of calculating the time derivative of the neutrino dis- 
tribution function is given in Scc.[n] The main result from Sec. [Q] is then used in Sec. 11111 to compute the neutrino 
emissivity. The specific heat is calculated in Sec. II VI Both quantities are used in Sec.[V]to discuss the cooling behavior 
of the considered phases. In Sec. IVII we explain the asymmetry of the neutrino emission in the A phase. Finally, we 
comment briefly on the effect of nonzero quark masses in Sec. I VIII 

Our convention for the metric tensor is = diag(l, — f , — f , — f ). Our units are h = c = Ub = 1- Four-vectors are 
denoted by capital letters, K = — (fco,k), while k = |k| and k = k/fc. We work in the imaginary-time formalism, 
i.e., T/VJ2k = TY^nf ^ 3 k/(27r) 3 , where n labels the Matsubara frequencies uj n = iko. The Matsubara frequencies 
are w„ = 2nirT for bosons, and w„ = (2n + l)irT for fcrmions. 

II. TIME DERIVATIVE OF THE NEUTRINO DISTRIBUTION FUNCTION 

In this section, we derive a general expression for the time derivative of the neutrino distribution function in spin-one 
color-superconducting phases. 

A. General formalism 

Within the Kadanoff-Baym formalism ppl |. one derives the following kinetic equation for the neutrino Green function 
(see for example Refs. HHMEl): 

id^Tr[ lx G<(X,P v )] = -Tr [G>(X,P v )i:<(X,P v ) - £> (X , P V )G< {X , P v )] , (1) 

where the trace runs over Dirac space. (Here and in the following, the index v always labels neutrino quantities and 
should not be confused with a Lorentz index.) The kinetic equation is obtained from the general Kadanoff-Baym 
equation after applying a gradient expansion, which is valid when the neutrino Green functions G^ y> {X, P v ) and 
self-energies E^ ,> (X, P v ) are slowly varying functions of the space-time coordinate X = (t, x). For our purposes, it is 
sufficient to consider spatially homogeneous systems which are close to equilibrium. In this case, the neutrino Green 
functions assume the following approximate form |22l |: 

-{-f X P v ,\ + A^7o)— {fu(t,p v )S{p v0 + [i v -p v ) - [I - fv(t, -Pu)] S(p v o + Hv+Pv)} , (2a) 

Pv 

{l X Pu,\ + P-vlo)— {[1 - fu(t, Pv)} S(puO + Vv - Pv) - fv(t, -Pv) S(PvO + fJ-v + Pv)} , (2b) 
Pv 

where /i„ is the neutrino chemical potential. The functions f V v{t, Pv) are the neutrino and antincutrino distribution 
functions. We are interested in the Urea processes u + e~ — > d+ v (electron capture) and d — > u + e~ + V (/3-decay), 
which provide the dominant cooling mechanism for quark matter. These processes are important for neutron stars 
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FIG. 1: Neutrino self-energies relevant for the neutrino Urea processes in the close-time-path formalism. The + and — signs 
assign the vertices to the upper and lower branch of the time contour, respectively. 



with core temperatures of the order of or smaller than several MeV. In this stage of the stellar evolution, the mean 
free path of neutrinos becomes larger than the stellar radius, wherefore in our results we shall set \i v = 0. 

The leading order contributions to the neutrino self-energies which enter the kinetic equation Q are given by the 
diagrams in Fig. ^ For the sake of simplicity, we do not take strange quarks into account. Their weak interactions 
are Cabibbo suppressed, and their number density is not expected to be very large. (Admittedly, however, the bigger 
phase space for the Urea processes involving massive strange quarks may partially compensate this suppression.) The 
diagrams in Fig. ^ translate into the following expressions, 

S<(*,P„) = J |^7 A (l-75)(7 K ^+Me7o)7 CT (l-75) 

X H> a (P e - P u ) — f e (t, Pe)S{ P eO + fl e - Pe), (3a) 
Pe 

Z>(t,P u ) = 7^7 A (l-75)(7 K ^e, K + Me7o)7 CT (l-75) 

X n< CT (P e - P v ) — [1 - f e (t,p e )]S( PeQ + M e -Pe), (3b) 
Pe 

where Gf is the Fermi coupling constant and fi e is the electron chemical potential. In the derivation we used the 
explicit form of the electron Green functions G^^it, P e ). They are given by expressions analogous to Eqs. We 
have neglected the positron contribution, i.e., the analogues of the second terms in curly brackets on the right-hand 
sides of Eqs. iJSJ. [Note that the processes involving positrons, u — > d + e + + v and d + e + — > u + 77, are suppressed 
by a large factor exp(2/z e /T) because /x e is positive and, in the regime under consideration, much larger than the 
temperature.] 

The functions n< CT > (P e - P v ) are the self-energies of the W bosons. The W exchange is approximated by its local 
form since the typical neutrino energies are much smaller than the W mass. In the case of neutrino processes in 
compact stars, the corresponding neutrino energies do not exceed several dozens MeV. 

Next, we insert the Green functions © and self-energies (0 into the kinetic equation (|TJ. In order to obtain the 
expression for the neutrino (antineutrino) distribution function, we integrate on both sides of the kinetic equation 
over p„o from — [i v to oo (from — oo to —[iv). The results are 

^/v(t,P„) = J {2 ^ p e vp L^(p e , Pv ) {[l-f v (t,Pu)]fe(t,Pe)n> a (Q) 



- f u (t, PlJ )[l- f e (t,p e )]U< a (Q)} , (4a) 

Mt, Pv ) = f ( /.f e L Ag (p e ,P,) {[l-/^, Pl/ )][l-/ e (i,Pe)]n< CT (Q') 

- fu{t,p v ) fe(t,Pe)^{Q')} . ( 4b ) 



with the four-momenta Q = (p e — p v — fi e + fi v , p e — p„) and Q' = (p e + p v — fi e + \i v , p e + p v ), and the shorthand 
notation 



L ACT (p e , p u ) = Tr [( loPe - 7 • Pe) 7 CT (1 - 7 5 )(7oP, - 7 • Vu) 7 A (1 - 7 5 )] 



(5) 
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phase 




M k 


Ak,i (m) 


Ak,2 (nz) 


Ak,3 (ns) 


CSL 


Sij 


J-7x(k) 


2(8) 


0(4) 




planar 


SilSji + Si2Sj2 


■/i7x,i( k ) + <^m,2(k) 


1 + cos 2 <9 k (8) 


0(4) 




polar 


5i3Sj3 


J37±,3(k) 


sin 2 9 k (8) 


0(4) 




A 


Sis(Sji + i 8j2) 


Js[7x,i(k) +«7x,2(k)] 


(l + |cos0 k |) 2 (4) 


(l-|cos0 k |) 2 (4) 


(4) 



TABLE I: Matrices A and A-^k and eigenvalues Ak,r with corresponding degeneracies n r for four spin-one color superconductors. 
The angle between k and the z-axis is denoted by 9 k . 



In the following, it is convenient to express the results in terms of the retarded self-energy for W bosons, LT^ 7 . 
Therefore, we shall use the following relations: 

n>(Q) = -2i[l + n B (a,)]ImIIfl(Q) , (6a) 
n<(Q) - ~2in B (q )ImTl R (Q) , (6b) 

where tib(lo) = l/(e w / T — 1) is the Bose-Einstein distribution function. 

Here we consider the cooling of quark matter only in the absence of neutrino trapping. Then, the neutrino and 
antincutrino distribution functions on the right-hand side of Eq. arc vanishing, f v .v(t, p„) = 0. The electron 
distribution function can be approximated by its equilibrium expression, 

fe(t, Pe) ^ n F (Pe ~ Me) , (7) 

with np{uj) = l/{e u l T + 1). Then, one arrives at 
d G 2 f d?p e 

-^S v {t,p v ) = 7n~yr~ — Lxa(p e ,Pu)n F (p e - He)n B {p v + Me - p e ) Imll^ (Q), (8a) 
at 8 J (ZirjJpvPe 

d G 2 f cPp 

■KlM^Pv) = r ro \* e — L x< y{Pe,Pu)n F {^ e -p e )n B (pu - Me + p e )lmlljf (Q')- ( 8b ) 

We shall see in Sec. Ill L>1 that the right-hand sides of Eqs. I|8a|) and (|8b|) arc in fact identical. 



B. Quark propagators in spin-one color superconductors 



As we shall see below, cf. Eq. I|18|l . the expression for the imaginary part of the polarization tensor of the W- 
vector boson is given in terms of the quark propagator S(K). In this paper, we consider the following four spin-one 
superconducting phases: the CSL, planar, polar, and A phases. Each of these phases is characterized by a 12 x 12 
gap matrix A'fk in color and Dirac space, 



(9) 



where (Ji)jk = —icijk and 7j_(k) = 7 — k^ k are basis vectors for the color antitriplet and spin triplet representations, 
respectively. We focus exclusively on the transverse phases, which have the highest pressure at asymptotic density 
and thus are expected to be preferred over others. For a more general form of the matrix A^k see Ref. ^?|- The 
matrix A is a complex 3x3 matrix and assumes a specific structure for each of the phases. In Table HJ we give A 
and the resulting matrices A^k- 

In spin-one color superconductors, the quark propagator is diagonal in flavor space, S{K) = diag[S u (K), Sd{K )). 
The Nambu-Gorkov structure of the flavor-diagonal elements is given by 



S f{ K)=( G l {K) ^ {K) 



f = u,d, 



where 01 [11 03 



Gf(K)= [G ^(A')] _1 



er ^0 ( £ k,r,/) 2 



(10) 



(11) 
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For the sake of simplicity, we have not included the quark self-energy correction. For a more general expression of the 
propagator, including this correction, see Refs. [2J,|2f|. The inverse free propagator for quarks and charge-conjugate 
quarks in the ultrarelativistic limit is 

G± f (K)] ^ = 7^ ± mo = 7o ]>>o ± (m/ - efc)] A± e . (12) 

e 

The Dirac matrices A£ = (1 + e7o7 ■ k)/2, where e = ±, are projectors onto positive and negative energy states. 
The matrices VZ r and V£ in Eq. are projectors onto the eigenspaces of the matrices Ai^M^. and 7°.Mj i .Mk7 t \ 
respectively. Both matrices have the same set of eigenvalues Ak r , 

M k Ml = 5> k , r P k - r , (13) 

r 

7 -M k A4 k7 ee ^A k>r P+ r . (14) 

r 

As we shall see, only the projection operators V^ r (and not V^ r ) are needed in the calculation of the neutrino emission. 
They are given explicitly in Appendix [BJ The eigenvalues A kr . appear in the quasiparticle dispersion relations, 



v /( e fc- M /) 2 + A k , r |0 / | 2 . (15) 



Here 4>f are the gap parameters, which are different for u and d quarks in general. The eigenvalues for the four 
considered phases as well as their degeneracies n r are listed in Table HJ Note that all phases contain an ungapped 
mode, A k , 2 (3) = 0. 

The off-diagonal elements on the right-hand side of Eq. (|10f> are the so-called anomalous propagators. They are 
given by 

E J( R ) = -E^^ k >f , V ■ ( 16b ) 

e,r k,rj) 

The propagators in Eq. (|11|) have particle and hole type poles at fc = ±e k r p as well as the corresponding antiparticle 
poles at fco = ±e^ r j. In the calculation of the imaginary part of the retarded self-energy H^(Q), the antiparticle 
contributions are suppressed by inverse powers of the quark chemical potential. [Note, however, that taking antipar- 
ticles into account is important in the calculation of ReH'^ T (Q).] In our calculation, therefore, we omit the terms with 
e = — in the quark propagator and arrive at the following approximate form 



^^^^ 7T >j , (17) 



where we denoted £k,r,/ = £ kr 



C. Imaginary part of the W-boson polarization tensor 

In this subsection, we evaluate Imll^ 7 . As can be seen from Fig.^ the W^-boson polarization tensor can be written 

as 

n ACT (Q) ee ^Tr[r^(^)r;5(A- + Q)] , (18) 

K 

where the trace runs over Nambu-Gorkov, color, flavor, and Dirac space. The 2x2 Nambu-Gorkov structure of the 
vertices reads 

I ^=Ml-7 B )r ± \ (19) 

± I - 7 a ( i +7 5 )tt J 
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where t± — (ri ± zt 2 )/2 are matrices in flavor space, constructed from the Pauli matrices Ti, r-x- After performing 
the traces over Nambu-Gorkov and flavor space, 

n Aff (Q) = ^£{Tr[ 7 \l- 7 5 )GM 7 CT (l- 7 5 )G+^^ , (20) 

K 

where P = K + Q, and the traces run over color and Dirac space. As we see, the anomalous propagators do not 
contribute to the self-energy. This is in contrast to spin-zero color superconductors, such as the (gapless) 2SC and 
CFL phases, in which the anomalous propagators are off-diagonal in flavor space. In fact, this difference is related 
to the conservation of electric charge in the diagrams in Fig. ^ The anomalous propagators contain the Cooper 
pair condensate, which, in the case of a spin-one color superconductor, is of the form (uu) or (dd). Therefore, the 
appearance of the anomalous propagators in the diagrams of Fig. ^ is forbidden by electric charge conservation. In a 
spin- zero color superconductor, however, the condensate is of the form (ud). Hence, electric charge can be extracted 
from or deposited into the condensate of Cooper pairs in the ground state and diagrams containing the anomalous 
propagators contribute to the W-boson polarization tensor [2(|. 

After inserting Eq. I|17|) into Eq. (|20|l . one arrives at the following expression for the imaginary part, 

TrnnWm - T T™ \^ \^ / t fc ° Z (gg Z k ^ Z ^ d ~ P ^ ft ftl , l k + (gj - bo + - P)] n-\a - A 

lmll fl - — im 2^2^) 7U "2 T7T2 72 \ J rs,+ { k ,P) + 7p ~2 wl2 ~2 A i r S ,-l K )PJ f . 

K rs { { £ k,r. J \Po ~ e P ,s.d) ( K ~ e U.r,d) \Po ~ e p,s,u) ) 

(21) 

where we defined the following traces in color and Dirac space, 

3# ± (k, p) = Tr [ 7 A (1 T 7 5 ) 7° 7>£ r AJ -f(l T 7 5 ) 7° ^ Aj] . (22) 
In order to perform the Matsubara sum, we use Eq. (|A4() in Appendix^ Then, extracting the imaginary part yields 

d 3 k 



t TT-Wm f d k \ x - »f(— eiek ,r ,«) n F (e 2 e p s d ) 

lmU R (Q) = -n^^J ^ [^(k.pl^Bp;,, ^(-e^ + e 2£p , s , g ^ ~ + e 2 e p , M ) 

(23) 



+ 7^f_(k,p) B k M B; ! Sj „ r - J <J(go + eie k ,r,d - e 2 e P)S ,„) 

Here, we defined the Bogoliubov coefficients 

B l,rj = ^ . {f = u,d, e = ±) . (24) 

The two terms on the right-hand side of Eq. i|23|) yield the same contribution. The physical reason for this is that the 
second term is just the charge-conjugate counterpart of the first term. The formal proof goes as follows. In the second 
term, one changes the summation indices ei <-> e 2 and r <-> s. Then, one introduces the new integration variable 
k — ► k — q and uses 

T s A r CT _(-p,-k)=T r ^ + (k,p) , (25) 

which holds for all phases considered in this paper. After taking into account that A k , 7 . = ^-k,r, one obtains the first 
term in Eq. (|23|) . Consequently, in the following we keep only the first term and double the result, 

d 3 k 



/(7 J k 
T^ + (k,p)B^ u B;^ d 
\ J 

S(q - eie k , r ,« + e 2 e P , s ,d) ■ (26) 



r,s ei ,62— i ' 

JIf(— eiek,r,u) n F (e 2 e p , B ,d) 



ns(— eiek,r,u + e2e p , s ,d) 

The same result may also be written in another form which is more convenient for the use in Eq. JsEJl, i.e 

/i7 3 k 
« 7# + (k,p')B£ PfU B 



r,s ei,e2=± 

n.F(eiek,r,«) np{-e 2 e p 



«s(eiek,r,« — e 2 e p ', s ,<i 
where we have used Eq. (|A2|) and defined P' = K + Q' . 



(27r) 3 rs >+ v ' r y k,r,« p',s,. 

5(9o - eie k ,r,« + e 2 e p ', s ,d) . (27) 
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D. General result for the time derivative of the neutrino distribution function 



We now insert the results for the polarization tensors (|26|l and Ij27[> into Eqs. (|8a|l and l|8b[) . respectively. In order to 
calculate L ACT (p e , p^) Imll^f^Q), we have to compute the contraction of the tensor L Ao .(p e ,p„) with the color-Dirac 
trace 7^°, (k, p). In all cases considered in this paper, we can write the result as 

WPe, P.) T r X s a +(k p) = 64 (p e - p e • k) {p v - p v • p) u r3 (k, p) , (28) 

where the functions aj rs (k, p) depend on the specific phase. They are calculated in Appendix 151 
Then, Eqs. lj%a[) and l|8b)l become 



-/,(t,p.) = -16nG% J J^y^j ^^(Pe-Me) £ £(Pe - Pe ■ k) (p„ - P„ - P ) 



ei,e2 — =b rs 



x u; rs (k,p) Bp^ d n f (-ei6 kir ^) n f (e 2 e P:S! ,i) <S(p e - {x e -p„- e^e^u + e 2 e p , s , £ z) , (29a) 

x w rs (k,p') ^ T . iU ^ iS d n F (eie kirjU ) <ii?(-e 2 ep, iM ) <5(p e - fi e + p u - eie kjn „ + e 2 e pVjd ) . (29b) 

At this point it is appropriate to recall the known result that the neutrino emission from low-temperature weakly 
interacting matter of massless quarks is strongly suppressed by the kinematics of the Urea processes [27J. In particular, 
energy and momentum conservation requires nearly collinear momenta of the participating electron, up quark and 
down quark (after taking into account that p v ~ T <C /U e , fx u , fid)- Strong interaction between quarks changes the 
situation dramatically [271 ]. In this case, applying Landau's theory of Fermi liquids, the quark Fermi velocity is 
reduced, vf — 1 — k, where k = 2a s /(37r) with the strong coupling constant a s . (Here we ignore non- Fermi liquid 
corrections. For their possible effects on the Urea processes in ungapped nuclear or quark matter see Refs. and pflj . 
respectively.) A rigorous treatment of the Fermi liquid correction would require a nonzero quark self-energy, which 
we have omitted, see Eq. <|11|) and remark below that equation. However, we may introduce the modified dispersion 
relations by hand in order to reproduce the result for normal quark matter as a limit case of our expressions. To 
this end, we identify the first three factors (multiplied by kp) on the right hand side of Eq. I|28|) with the squared 
scattering amplitude |M| 2 , derived by Iwamoto [271 ]. 

64G F ( Pe k - Pe • k)(p„p - p v ■ p) -> |M| 2 = MG 2 F (P e • K){P V ■ P) . (30) 

After this replacement, we take the Fermi liquid corrections into account by simply following the same steps as in 
Ref. [23] ■ As in Ref. [23, we work to lowest order in a s . Therefore, the results of this paper are valid, strictly speaking, 
only at densities much higher than in the interior of neutron stars. However, the role of the strong interaction in 
the neutrino emission is merely to open a phase space for the weak processes. In view of this, the limitation due to 
uncontrollable strong interaction may not be so essential for understanding the qualitative features of the neutrino 
processes in quark matter at realistic densities. 

We approximate the 6- functions in Eqs. (|29[l as follows. By making use of the definitions p = k + p e — p„ and 
p' ee k + p e + p„, it is easy to show that the arguments of the 6- functions vanish only if the angle between the 
momenta of up and down quarks 9 u d is equal to a fixed value 9q, up to corrections suppressed by powers of p, y //i e . 
The value of the angle 9 is given by cos^o = 1 — K/J%/(Li u [J,d)- We note that 9q is independent of the neutrino energy 
p v . Then, both 5-functions in Eqs. I[29a[) and l[29b[) are replaced by fi e / (^uHd) S(cos6 U( i — cos#o). Making use of this 
fact, we rewrite the equations in the following form, 

-g^Jiy(t,Pu) ^ -Q s G f J J ~(2^jT (1 ~ cos ^rf) 0(COs9 ud - COS 6> ) 

x X! w rs(k, p) B^ r u B p 2 s d n F (-eie ktr>u ) n F (e 2 e v , s ,d) n F (Pu + eie k>f . )U - e 2 e p , S: d) , (31a) 

ei,e2 — ± rs 

d f u \ 64 n-2 f d pdil p f dkdilu 

t^M^PW — — -a s Ci F fie^u/id J J (27r) 3 ^ ~ C ° S ' ^ C0S ~ C0S ' 

x ^u rs (k,p)B^ ru B^ sd n F (e 1 € Kr ^)nF{-e 2 e v , s ^)nF(p v ~ , (31b) 

ei,e2 — ± rs 
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where 9„d is the angle between the three-momenta of the neutrino and the d quark. Here we have changed the 
integration variable p e to p = k + p e — p„ and p' = k + p e + p u , respectively, and afterwards dropped the prime of 
p' in Eq. (j31b)) . 

Instead of dimensionful momenta, it is convenient to introduce the following three dimensionless variables, 

P-Vd fc-Mu Vv 

x = T , y= T , V = Y' ^ ' 

In terms of the variables x and y, the range of integration runs from —(i u ,d/T to oo. Since the main contribution in 
the integral comes from x,y fJ- u ,d/T, the results do not change if the lower boundary is extended to — oo. Then, 
we can drop the x- and y-odd contributions from the Bogoliubov coefficients, 

K's.d = \ - 7 2 6 f x 5 , (33a) 



K)r,u = \ Z , = i . 



in the integrand, and keep only the even contributions, i.e., a constant term 1/2. After this is taken into account, 
we restrict the integration over x and y from to oo. Accounting for the integration from — oo to produces extra 
factors of 2 for each integration which compensate the factors of 1/2 from the two Bogoliubov coefficients. At this 
point, we notice that the difference between Eqs. (|31a[) and (|31b() lies only in the signs in front of e\ and e^- Because 
of the summation over e\ and e^-, the right-hand sides of both Eqs. I|31a|) and (|31b() arc identical. 
Hence, the result can be written in the following approximate form, 

- -y a s G 2 F n E n u ti d T 2 Y^ J J (1 - cos6W) S(cos9 ud - cos6> ) F^ d {% p, v) , (34) 

where 



F£ Vd fcp,v) = o; rs (k, P ) J] / / dxdy(e- e ^y 2 +^l + i) ( e *\/^+W5+l) 

x ^ e f+ei-v/« 2 + A k.'-vi-e2- v / a;2 + A P.= ¥ , 3 + 1^ 1 (35) 

with 93/ = 4>f/T. 

Three out of four angular integrations in Eq. (|34|l can be calculated approximately in an analytical form. Details 
of this calculation are deferred to Appendix IO We obtain the following main result of this section, 

Mt,P») = §- t Mt,P»)--^^^uVdT 2 Y, J dt(l-ZcaBe v )F£ Vd (Z,v) . (36) 



d_ 

di 



It holds for isotropic phases as well as for phases in which the order parameter picks a special direction in momentum 
space, identified with the z direction. We denote the angle between the neutrino momentum and the z-axis by B v . 
The function F£,* (j£, v) is obtained from the function i^™ (k, p, v) in the collinear limit (k = p). Its explicit form 
reads 



/>oo />oo _1 . . _i 

*Z Vd &v) = WrrK) E / / dxdy(e-^V^h^l+l) (W* 2 +W,+i) 

e u e 2 = ± J JO V ^ V ' 

Here £ = cos U = cos 6^ with 8 U being the angle between the z-axis and the u quark momentum k, and 8d being the 
angle between the z-axis and the d quark momentum p. The functions w rr (£) and A^. r are given in Table ITT1 One 
should note that all components of o; rs (k, p) with r 7^ s disappear in the collinear limit. This can be seen directly 
from their explicit expressions given in Appendix [5] Because of this property, every excitation branch r yields a 
separate contribution to the result 136|) . 
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phase 






1^22 (£) 


A{,2 


^33(C) 


Aj,3 


CSL 


2 


2 


1 









planar 


2 


l + C 2 


1 









polar 


2 




1 









A 


1 + sgn(C) 


(1 + lel) 2 


1 - sgn(C) 


(i-ICI) 2 


1 






TABLE II: Functions cj rr (£) and A^ >r for four spin-one color superconductors. 



III. NEUTRINO EMISSIVITY 



In this section, we calculate the neutrino emissivity. 

d 3 



(38) 



This is the total energy loss per unit time and unit volume carried away from quark matter by neutrinos and 
antineutrinos. Inserting Eq. (|36[) into Eq. I|38|) and making use of the integral 



E 



dv v" 



dx 



dy (e- eiv + 1) 1 (e e2X + l)" 1 [ e v+eiV ~ e2X 



' 5040 



(39) 



we obtain 



where 



630 



1 2 

2 + ^G(p u ,p d ) 



G{Vu,Vd) 



1260 

457tt 6 j 
1260 



457tt 6 



dvv 3 
dv v 3 



(CSL, planar, polar) . 



(40) 

(41a) 
(41b) 



In all phases we consider, the emissivity t v consists of two contributions. The first contribution is given by the term 
1/3 in the square brackets on the right-hand side of Eq. Q40[l. It originates from the ungapped modes: r = 2 in the 
CSL, planar, and polar phases, and r = 3 in the A phase. The second contribution is given by the term proportional 
to G(cp u , ifd)- It originates from the gapped modes. The function G(ip u , ifid) has to be evaluated numerically for each 
phase separately. For the sake of simplicity, we set tp u = (fd = (p in the following. The results for G(tp, <p) are plotted 
in the left panel of Fig. |21 The right panel of Fig. [21 shows the function G(T), obtained from G(<p, <p) by making use 
of the following model temperature dependence of the gap parameter, 



( y 



(42) 



with 0o being the value of the gap parameter at T = 0, and T c being the value of the critical temperature. 

As a consistency check, we first read off from the figure that the general result in Eqs. (|4(J|) and (|41|) reproduces the 
well-known expression for the neutrino emissivity in the normal phase. This is obtained by taking the limit ip —> 0. 
Of course, the result in this limit is the same for all considered phases. Since G(0,0) = 1, see left panel of the figure, 
we recover Iwamoto's result p7T |. 

In the spin-one phases, the function G((p, ip) describes the suppression of the emissivity due to the presence of the 
gap in the quasiparticle spectrum. In order to discuss this suppression for the different phases, we derive analytic 
results for the asymptotic behavior of G(p, ip) for p — > oo. Physically, this corresponds to the low temperature 
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FIG. 2: Left panel: the suppression functions G(f>, <p) of the neutrino emission contributions due to gapped modes in the CSL, 
planar, polar and A phases. Right panel: the temperature dependence of the suppression functions using Eq. 1421 . 



G(<p,<p) ~ < 



1 in Appendix |D| The 




(CSL) , 




(planar) , 




(polar) , 




(A). 



(43) 



The strongest suppression happens in the CSL phase, in which the gap is isotropic. At large values of ip, the 
cmissivity is exponentially suppressed, which is universal and qualitatively the same as, for example, in a spin-zero 
color superconductor [26J or in superfluid neutron and/or proton matter [30l l3lj. From a physical viewpoint, this 
reflects the fact that, in a superconductor, the neutrino emission is proportional to the density of thermally excited 
quasiparticles. It is worth emphasizing, however, that the function G(ip, <p) cannot be approximated well by the 
exponential function at small (p. For tp < 1, the actual suppression is much weaker. This is also obvious from the 
right panel of Fig. [3 This panel shows that, in all phases, the function G(T/T C ) behaves almost linearly all the way 
down to T/T c ~ 0.4. In the CSL and the planar phases, the exponential suppression starts to show up only below 
this point. 

The contributions from the gapped modes in the other spin-1 phases differ considerably from the CSL result. All 
of them have some degree of anisotropy in the gap function. The second strongest suppression is seen in the planar 
phase. In this case, while the gap function is anisotropic, it has no zeros. The dominant contribution comes from a 
stripe around the equator of the Fermi sphere, where the gap function, which is proportional to */! + cos 2 / , takes 
its minimum. 

In the polar phase, the gap function has point nodes at the north and south poles of the Fermi sphere, i.e., it costs 
no energy to excite quasiparticles around these points. Therefore, these quasiparticles give the dominant contribution 
to the emissivity. For large values of the dimcnsionlcss variable <p, in particular, one has a power-law instead of the 
exponential suppression. 

The gap function in the A phase has also nodes at the north and south poles of the Fermi sphere. However, there 
is a difference compared to the polar phase in the behavior of the dispersion relations at small angles / . While it is 
linear in Of in the polar phase, it is quadratic in Of in the A phase. This difference gives rise to a different suppression 
of the emissivity, see Fig. |3 and Eq. I|43|) . 

The results of this section will be used in Sec. in order to discuss the effect of spin-1 superconductivity on the 
cooling rates of compact stars. Besides the neutrino cmissivity, this requires the calculation of the specific heat. 
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IV. SPECIFIC HEAT 



In this section, we calculate the specific heat of spin-one color superconductors. This result shall be used in the 
discussion of the cooling rate in the next section. We may start from the entropy density (see for example Ref. |l8j| ) 



/73 
T^jTg {nF(ek,r,/) lnnF(ek,r,/) + [1 - n F (e KrJ )\ ln[l - n F (ek,r,/)]} , 



(44) 



f=u,d r 



where n r is the degeneracy of the quasiparticle branch r as given in Tabled I n each phase, ^ r n r /2 = 6, accounting 
for 2 spin and 3 color degrees of freedom. The specific heat is then obtained as 



T 



dS_ 
df 



EE 

f=u,d r 



d 3 k d 
—^£k,rj-g^n F (e KrJ ) 



(2tt) s 



(45) 



Making use of the model temperature dependence for the gap parameter in Eq. (|4*2*|) . the result for all phases we 
consider can be written as 



c v = T E A 

f=u,d 



(46) 



The structure of this result is analogous to that of emissivity in Eq. (|40|) . i.e., the first and the second terms in the 
square brackets on the right-hand side come from ungapped and gapped modes, respectively The explicit form of the 
function K(ip) reads 



K(<p) 



IT 2 , 



dx 



d Z 9 

i L^+x^r* + ^ 



x z + A e ,i <p 



dx 



di 



X Z + Af.r ( V + 



(CSL, planar, polar), (47a) 
(A) . (47b) 



The numerical results for the function K(ip) for all considered cases are plotted in the left panel of Fig. [3J The 
value of this function in the limit ip — > does not reproduce the result for the normal phase K n = 1, we rather 
observe lim^^o K(ip) > 1. This is due to the jump of the specific heat Acy at the point of the second order phase 
transition to superconducting matter. In the model used, the magnitude of the jump is obtained by calculating the 
term proportional to </>§/T c 2 in Eqs. (|47|l . The result can be written as 



,-r2 



Ac v 
Ac v 



{A + pi) 



(CSL, planar, polar) 
(A) , 



(48a) 
(48b) 



where 4> r = (f dCl^/ (Att) Ak.r) 1 ^ 2 is the quadratic mean of the r-th gap function. Using the results from Ref. [l7j . 
cf. Eq. (118) therein, we conclude that the jump of the specific heat is proportional to the condensation energy (at 
T = 0). Hence the order of the values K (0) in Fig. |3| reflects the order of the condensation energies. 

As for the emissivities, we derive analytical approximate expressions for the specific heat at asymptotically large 
ip, corresponding to asymptotically small temperatures. For the details of the calculation see Appendix El We find 



K(<p) 



^5/2 e -V2 V 


(CSL) , 


p 2 e-* 


(planar) 




(polar) , 


p- 1 


(A). 



(49) 
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FIG. 3: The functions i£"(y>) (left panel) and K(T/T C ) (right panel) for four spin-one color superconductors. 

This behavior of the specific heat leads to the fact that the curves in Fig. reverse their order at large <p compared 
to ip = 0, e.g., the specific heat in the CSL phase, which has the largest jump at the critical temperature, exhibits the 
largest suppression for very small temperatures. 

Note that the degree of the suppression due to gapped modes in the specific heat and in the emissivity are similar 
at large tp, cf. Eqs. and 

V. COOLING RATES 



In this section, we shall use the results for the emissivity in Eq. (|40Jl the specific heat in Eq. (|46() in order to 
study cooling of bulk matter in spin-one color-superconducting phases. When the cooling is only due to the neutrino 
emissivity, one has the following relation, 



e„(T) = -c v (T) 



dT 
~dt 



In order to derive the change of temperature in time, one has to integrate the above equation, 



t-tn 



T dT> C -^l 
t e„(T') ' 



(50) 



(51) 



where Tq is the temperature at time to. By inserting the expressions from Eqs. (|40() and i|46|) into Eq. (|51|l and using 
<Pu — = ¥>! we derive 



t-tn = 



630 



/4 



dT' 1 + 2K(T') 
457 a s G 2 F n e fi u n d J ' To (T') 5 1 + 2G(T") 



(52) 



where the temperature-dependent functions K(T) and G(T) are obtained from the functions K{ip) and G(ip, p) with 
the help of Eq. 

By making use of Eq. I|52|l . let us estimate the cooling behavior of a compact star whose core is made out of spin- 
one color-superconducting quark matter. We start from the moment when the stellar core, to a good approximation, 
becomes isothermal. At this point, the stellar a ge i s of the order of £ = 10 2 yr and the temperature is about 
To = 100 keV. The estimates in the literature |15l l32| suggest that the value of the critical temperature in spin-one 
color superconductors is of the order of T c — 50 keV. This is the value that we use in the numerical analysis. Moreover, 
we choose /x u = 400 MeV, \i& = 500 MeV, /x e = 100 MeV, a s = 1. The Fermi weak coupling constant is given by 
G F = 1.16637 x 10" 11 MeV" 2 . 

The numerical results show that the cooling behavior is dominated by the ungapped modes. Consequently, to a very 
good approximation, the time dependence of the temperature can be computed by neglecting the functions K(T') 
and G(T') in Eq. (|52|) . In this case, an analytical expression can be easily derived, 



T(t) 



Tpr 1 / 4 
(i-i + r)i/4 ' 



(53) 
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FIG. 4: Temperature as a function of time for normal quark matter and four spin-one "toy phases" (dropping the ungapped 
modes). The curves represent the CSL phase (solid), planar phase (long-dashed), polar phase (short-dashed), A phase (dotted), 
and normal quark matter (dashed-dotted) . 



where 

r = 315 &±ti 1 (54) 
914 a s G% ^ (J-u Tq 

With the above parameters, this constant is of the order of several minutes, r ~ 10~ 5 yr. 

It may be interesting, although unphysical, to compare the cooling behavior of the gapped modes of the different 
spin-one phases. To this end, we drop the 1 in the numerator and denominator of the integrand in Eq. Ij52[l. The 
results are shown in Fig.^J Note that both the initial temperature T and the critical temperature T c are beyond the 
scale of the figure. The reason is that, even for the gapped modes, the cooling time scale for temperatures down to 
approximately 10 keV is set by the above constant r. Therefore, all phases cool down very fast, and the transition to 
the superconducting phase at T — 50 keV is hidden in the almost vertical shape of the curve. Only at temperatures 
several times smaller than T c , i.e., of the order of 10 keV, substantial differences between the phases appear. In this 
range, the fully gapped phases cool down considerably slower than the phases with nodes on the Fermi sphere, which, 
in turn, cool slower than the normal phase. It seems to agree with physical intuition that this order reflects the order 
of the suppression at low temperatures for the neutrino emissivity, i.e., the slowest cooling (isotropic gap) happens 
for the phase where e v is suppressed strongest while the fastest cooling (no gap) happens for the smallest suppression. 
Note, however, that the cooling depends on the ratio of the suppressions of e„ and cy. Therefore, this order is a 
nontrivial consequence of the exact forms of the functions G(tp, tp) and K(tp). For large values of (p, we may use Eqs. 
(|43|1 and l|49|l to estimate the ratio K(ip)/G(ip, <p). For both completely gapped phases we find K(ip)/G(tp, <p) ~ tp 3 / 2 
while both phases with point nodes yield ratios independent of (p. Consequently, for late times, T ~ t~ 2 l xl in the 
CSL and planar phases while T ~ t~ x l A in the polar, A and normal phases. 



VI. SPATIAL ASYMMETRY IN THE NEUTRINO EMISSION FROM THE A PHASE 



In this section, we address a special aspect of the angular distribution of the neutrino emission. To this end, we 
consider the net momentum carried away by neutrinos and antincutrinos from the quark system per unit volume and 
time, 

rfP ( " et) d [ d 3 p„ d f d 3 p v 

-dvdf ^-dijj^r^ [/ ^' p - } + Mt > p » )] = ~ W w ^ ^ Vv) ■ (55) 

Analogously to the case of the total emissivity, see Eq. (|40|) . we arrive at the following general result, 

d p(net) 457 



dV dt 945 



a s G F T He^uiXd H((p u , <p d ) z , (56) 
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FIG. 5: Numerical results for the functions H(p,tp) and H(T/T C ) which determine the net momentum carried away from the 
spin-1 color superconducting A phase by neutrinos. 



where z is the unit vector in z direction, and 



420 f°° f 1 

H(tp u ,tp d ) = —^g J dvv 3 J dttF£ Vi (£,v) = (CSL, planar, polar) , (57a) 



420 r°° f 

H(tp u ,tp d ) = J dvv 3 J dZZ[F£ Vd (Z,v)+F? uVi (Z,vj\ (A). (57b) 

In the CSL, planar, and polar phases, the function H(tp Ul tp d ) is identically zero. This is because F^ Vd (t;,v) is an 
even function of £ in these three phases, and therefore the integration over £ in Eq. Ij57a(l is vanishing. This means 
that the net momentum of emitted neutrinos as well as the related net recoil momentum of bulk quark matter in the 
CSL, planar, and polar phases are zero. 

The result is non-vanishing, however, in the A phase. The corresponding function H(tp, tp) is plotted in the left 
panel of Fig. [5] From the figure, we see that H(0, 0) = 0. Of course, this is just a consistency check that, in the limit 
tp — > 0, we reproduce the vanishing result in the fully isotropic normal phase of quark matter. From the numerical 
data, we find that the maximum value of the function H (tp, tp) is approximately equal to 0.064, which corresponds to 
the value of its argument tp ~ 2.9. At large tp, the asymptotic behavior of H(tp,tp) is power suppressed as 1/tp. For 
completeness, we also show the temperature-dependent function H(T) in the right panel of Fig. [S] This is obtained 
from H(tp, tp) by making use of the model temperature dependence for the gap in Eq. 1421) . 

It may look surprising that the net momentum from the A phase is nonzero, indicating an asymmetry in the 
neutrino emission with respect to the reflection of the z-axis. The gap functions do not exhibit this asymmetry. 

The origin of this remarkable result can be made transparent by rewriting the expression (|57b|) for the A phase in 
the following way, 



840 r°° 

= ~ J o dvv 3 j d£t Ff uVd (£, v) , (58) 



where 



OO pOO 



x ^ e v+eWy 2 + (l+i) 2 vl-e2y/x 2 + (l+i) 2 <P 2 d + lj^ . (59) 

In the derivation, we used the explicit forms of ov r (£) and X^. r from Table ITU Now, the result looks as if only one 
single quasiparticle mode contributes to the net neutrino momentum. The corresponding "effective" gap function has 
the angular dependence ~ (1 +£) which clearly discriminates between +z and —z directions. 

In order to understand the physical reason for the appearance of the effective quasiparticle mode, it is useful to 
analyze the physical properties of the gapped modes of the A phase, r = 1,2. The color-spin structure of these modes 
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FIG. 6: Gap functions for the first (left) and the second (middle) excitation branch with specified helicities of quasiparticles 
in the upper and the lower hemispheres. The "effective" gap relevant for the neutrino emission is shown on the right. 

is encoded in the projection operators V£ . Their explicit form is given in Eq. (|B12fl . It is instructive to write the 
first two projectors in the form 

K.i = l J H l - s S n (^)] H+(k) + \ Ji[l + sgn(fc 3 )] F-(k) , (60a) 
V+ 2 = \ J 3 2 [l + sgn(fc 3 )] H+(k) + I J|[l - sgn(fc 3 )] H~(k) , (60b) 

where ff ± (k) = i (1 ± X • k) are the helicity projectors with X = 7 5 7°7. From Eq. I|60a|) we see that the quasiparticles 
of the first branch have helicity +1 when the projection of their momentum onto the z-axis is negative, fc 3 < 0, and 
helicity —1 if fc 3 > 0. Quasiparticles of the second branch have opposite helicities, see Eq. (|60bll . 

The next step in the argument is to notice that only left-handed quarks participate in the weak interactions which 
underly the Urea processes. Formally, this can be seen from Eq. I|20|l where the left chirality projectors |(1 — 7 5 ) occur 
in the first term under the trace. (Note that the second term describes charge-conjugate quarks for which i(l + 7 s ) 
projects also onto left chirality states.) In the ultrarelativistic limit, these are quarks with negative helicity. Taking 
into account the helicity properties of the quasiparticles in the A phase, it becomes clear that only an effective gap 
structure contributes. This is constructed from the upper hemisphere of the first mode and the lower hemisphere of 
the second mode, see Fig.El This is a graphical representation of the formal argument given after Eq. (|59|) . [Of course, 
our choice for the angular dependence of the gap functions, namely Ak,i = (1 + | cos#k|) 2 and Ak,2 = (1 — | cos#k|) 2 , 
is only one possible convention. Equivalently, one could choose Ak.i = (1 + cos#k) 2 and Ak,2 = (1 — cos6>k) 2 , in 
which case the quasiparticle excitations would be ordered according to their helicity. Then, quasiparticles of the first 
(second) branch would have negative (positive) helicity, and the weak interaction would involve only quasiparticles of 
the first branch. Our convention in this paper is in accordance with Ref. [TtI].] 

The asymmetry in the effective gap function translates into the asymmetry of the neutrino emission. This is due to 
the angular dependence of the amplitude for Urea type processes. As in the vacuum, the corresponding amplitude is 
proportional to 1 — cos9 u( i- This can be seen, for instance, from the integrand in Eq. I|34|) . Such an angular dependence 
of the amplitude means that the neutrinos are emitted preferably in the direction opposite to the (almost collinear) 
momenta of the participating up and down quarks. In fact, this is a general property that holds also in the normal 
phase |27j . Since the effective gap function assumes smaller values for quasiparticles with fc 3 < than with fc 3 > 0, 
there is more neutrino emission in the +z direction. 

One can estimate the maximum velocity of a neutron star with a quark matter core in the A phase that can be 
obtained by the asymmetric neutrino emission. It has been shown that this velocity is negligibly small, e.g., of the 
order 1 m/s, see erratum in Ref. ^ij^. In essence, the reason for this is that the available thermal energy in the star, 
after matter in the stellar interior cools down to the critical temperature T c < 100 keV of the A phase, is too small to 
power substantial momentum kicks. (It would be interesting to investigate, however, if additional sources of stellar 
heating, e.g., such as the latent heat from a first-order phase transition, could change the conclusion.) 

VII. QUARK MASS EFFECTS 

In nature, quarks are not exactly masslcss. Therefore, it is important to address the effects that the masses have on 
the dispersion relations of quasiparticles, and thus on the neutrino emissivity and the specific heat of quark matter. 
In order to study the massive case, we keep the color-spin structure of the gap matrix Alk exactly as in the massless 
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case, given in Eq. (j5J. Then, the inverse full quark propagator can be written as follows 

(61) 



with [Gj^(-ftT)] = l^K^ i 7°/i + to, where to is the quark mass. For simplicity, we omit the flavor index / in this 
section. 

We do not repeat the detailed calculations for the emissivity and the specific heat with this propagator. Instead 
we assume that the main modification happens due to the change of the quasiparticle dispersion relations. These 
relations are determined by the solution to the algebraic equation det S^ 1 = 0, or more explicitly, 



(let 



(G o -)- 1 (G+)- 1 -0 2 (G o )- 1 X k G o - 7 o X] c 7 = 0. (62) 



We find that, in all considered spin-one phases, the results for the dispersion relations are essentially the same as in 
the massless case, except for the replacement k — ► V k 2 + m 2 , i.e., 



e 



, 2 

2 
k. 



Vfc 2 + m 2 - ii ) + \ Kr <t> 2 , (63) 



with the same Ak,r's as for m = 0, sec Table [I] This result suggests that the emissivity and the specific heat do 
not change qualitatively after including quark masses. This conclusion may not be so surprising because, in general, 
non-zero quark masses are not expected to affect much the physical properties which are dominated by quasiparticle 
states in the vicinity of the Fermi sphere. 

In view of the above "trivial" effect of the quark masses, it is appropriate to comment on the recent study in 
Ref. j^il where it is argued that there are no ungapped modes in the CSL phase when the quarks are massive. This 
may look as a contradiction to our result showing that all ungapped modes, Ak,2 = 0, survive after switching on 
the mass. The seeming contradiction is removed, however, after noticing that a different choice of the gap matrix in 
the CSL phase is utilized in Ref. |34l |. In our notation, the corresponding gap matrix would be obtained by replacing 
7 i i (k) with 7j in Eq. (j5J. After making such a replacement, we find that the dispersion relations of quasiparticles in 
Ref. [24j are indeed reproduced. In particular, the low-energy dispersion relations for two out of total three different 
quasiparticles are given by 



,2 

t k,l/2 




(^W-„) ± " . (64) 



Note that these are obtained without the limitation of the smallness of the quark mass to. The corresponding two 
energy gaps are thus 




i/ 2 = V20 Wl+ — ±-_ . (65) 



The low-energy approximation, defined by £k,i/2// i ~ 4>l H ~ Wk 2 + m 2 — fi\/ fi <C 1, is completely sufficient for the 
study of most transport and neutrino processes in spin-one color-superconducting phases. 
The dispersion relation of the third quasiparticle mode can be extracted exactly, 



4,3 



/i 2 + k 2 + m 2 + <f) 2 — 1\/ fi 2 (k 2 + m 2 ) + k 2 <p 2 

m 2 4> 2 
~ 2 



(yk 2 + m 2 -fiy + ^-. (66) 



Hence the value of the energy gap is 



m(j) mcj) 



(67) 



For m = 0, we recover the gaps <j)\ = <p2 — V2 4>, <p3 = 0. Thus, in contrast to the massive case, both CSL gap matrices 
(i.e., one with 7j_j(k) and the other with jj in the definition of A4u) give rise to the same dispersion relations in the 
ultrarelativistic case. It should be studied in the future which of the two physically different CSL phases at to 7^ 
has the lower free energy. 
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VIII. CONCLUSIONS 

In this paper, we have computed the neutrino emissivity due to direct Urea processes (i.e., u + e~ — ► d + v and 
d — > it + e~ + p), as well as the specific heat in four different spin-one color-superconducting phases of dense quark 
matter. Starting from the kinetic equation, we have derived a general expression for the neutrino emissivity. In the 
case of the normal phase of quark matter, this reduces to the well-known analytical result j23|. The basic ingredients 
in the calculation are the quasiparticle dispersion relations, containing the spin-one gap functions. We have studied 
in detail the effect of an isotropic gap function (CSL phase) as well as of anisotropic gap functions (planar, polar, A 
phases). The numerical results for the emissivity and the specific heat as functions of the ratio of the gap parameter 
to the temperature, tp = <p/T, have been presented. 

In all four phases, also analytical expressions have been derived in the large tp limit (i.e., in the limit of small 
temperatures). In particular, in the case of an isotropic gap function (CSL phase), the well-known exponential 
suppression of the emissivity and the specific heat is observed. We find that anisotropic gaps give rise to different 
asymptotes in general. For example, the phases in which the gap has point nodes (i.e., polar and A phases) show 
a power-law instead of an exponential suppression at tp — ► oo. The actual form of the power-law depends on the 
behavior of the gap function in the vicinity of the nodes. While a linear behavior gives rise to a suppression ~ 1/y 2 , 
a quadratic behavior leads to ~ I/93. 

We have used the results for the emissivity and the specific heat to discuss the cooling curves of spin-one color 
superconductors. Our simplified analysis, performed for an infinite homogeneous system, reveals several important 
qualitative features. Most importantly, the cooling rates of all considered spin-one color superconducting phases differ 
very little from that of the normal phase. The reason is that, with the choice of the gap matrix as in Eq. |JSJ, all phases 
(even if the quarks are massive) have an ungapped quasiparticle mode which dominates cooling at low temperatures 
and, thus, makes the suppression effect due to gapped modes unobservable. Note, however, that this conclusion may 
change if another choice of the gap matrix, e.g., such as used in Ref. [34|], corresponds to the true ground state of 
quark matter. 

In addition to the cooling rates, we have discussed an unusual property of the color-superconducting A phase, in 
which the neutrino emission is not symmetric under the reflection of one of the coordinate axes in position space. This 
is seen from the fact that the net momentum of emitted neutrinos is nonzero, pointing into a direction spontaneously 
picked by the order parameter. As we have argued, the asymmetry is related to the helicity properties of the 
quasiparticlcs in the A phase. A helicity order arises naturally from the structure of the gap matrix which is a 
straightforward generalization of the A phase in 3 He (where, however, there is no helicity order). 

So far, we did not find any observable consequence of the helicity order in the A phase. The simplest possibility 
would be realized if the asymmetric neutrino emission could result in a "neutrino rocket" mechanism for stars 33]. 
The estimated effect on the stellar velocity appears to be extremely small, however. The search for other observable 
signatures, e.g., dealing with the timing of pulsars, may reveal other possibilities. We could also imagine that observ- 
able signatures of the helicity-typc order might be observed in completely different systems in atomic or condensed 
matter physics (e.g., trapped gases of cold fcrmionic atoms, or high-T c superconductors) if they happen to have a 
similar structure of the order parameter, see for example Ref. |35l |. A systematic study of such a possibility is outside 
the scope of this paper, however. 
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APPENDIX A: MATSUBARA SUM 



In the calculation of the polarization tensor Il : ^ 7 (Q) in Sec.|n]we use the following Matsubara sum, 



T E 



(fco+6)(fco + go+6) 
^(fc 2 -6?)[(fc + go) 2 -ei] 



4eie 2 



(ei- £0(62 + 6) (ei + £i)(e 2 -6) 



g - ei - (2 



go + ei + £2 



[1 - n F (ei 



n F (e 2 



(ei+£i)(e 2 +6) (ex-6)(ca-6) 



[n F (ei) - n F (e 2 )] 



go + ei - e 2 go - ei + £2 

where /co (go) is the fermionic (bosonic) energy and £i,£ 2 ,ei,e2 are real numbers (ei,£2 > 0). With 

n F (e 1 )n F (e 2 ) n F (-ei) n F (-e 2 ) 



and 



1 - n F {ei) - n F (e 2 ) 



n F (e 1 ) - n F (e 2 ) 



«s(ei+£2) ns(-ei-e 2 ) 
nF(-ei)nj?(e 2 ) n F (ei) n F (-e 2 ) 



ns(-ei+£2) n B (ex-£2) 
we can write the result in the following compact form, 

(fco + £i)(fco + go +6) 1 \ - (ei - ei^i) (e 2 - e 2 £ 2 ) n F (-eiei) n F (e 2 e 2 ) 



T E 



(fc 2 -e 2 )[(fc +g )2- 



E 



4eie 2 go-eiei+e 2 e 2 n B (-eiCi + e 2 e 2 ) 

ei,e 2 =± 



This formula is applied to Eq. l(2"T|) . 



(Al) 



(A2) 



(A3) 



(A4) 



APPENDIX B: COLOR AND DIRAC TRACES 

In this appendix, we compute the color and Dirac traces of the tensor 7^°, (k, p), see Eq. (12211 . in the polar, planar, 
A, and CSL phases. Moreover, we contract this tensor with the tensor L A<T (p e , p), see Eq. JjJ, to obtain the functions 
u! rs (k, p) , see Eq. (|28|) . In this appendix, for the sake of clarity, we shall add the subscripts c or D at the symbol 
"Tr" to indicate whether the trace is taken over color or Dirac space. In the calculation, we encounter the following 
expressions, 

-00/ 



L (Pe,Pu) = 8 (j> e Pv + p e • p„) , 

L 0i (p e ,p v ) = 8 \p e p % v +plPu + i (Pe X p v )i] , 

L M (p e ,Pv) = 8 \p e pl +plPu - i(Pe X P„)i] , 

£ ij '(p e ,P„) = 8 [5 ij (p ePv - p e ■ p„) + plpi + p 3 e pl -ie iji (p e pl - pip*)} 



where i, j = 1,2, 3, and 



r 00 (k, P ) = 2(i + k- P ) , 

T 0l (k,p) = 2[F+ ! i' + l (kxp)] , 

r i0 (k,p) = 2[k i +p i -i(kxp)} , 

T y '(k,p) = 2[5 ij (l-k-p) + k i p j + tip l + ie lji (k £ -f) 



where we abbreviated 



T ACT (k, p) = Tr D [ 7 A (1 - 7 5 ) 7 ° A k 7 CT (1 - 7 5 ) 7° A p ] 
Using the above results, we calculate the Lorentz contraction 

L\a{Pe,p v )T x,T (k,p) = 64 (p e - p e ■ k) (p v - p v p) 



(Bla) 
(Bib) 
(Blc) 
(Bid) 



(B2a) 
(B2b) 
(B2c) 
(B2d) 

(B3) 

(B4) 



This is a generic result which we shall use in the following subsections to compute the more complicated contractions 
for the spin-one phases. 
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1. Polar phase 



The polar phase is particularly simple, since the projection operators V^ r do not depend on the quark momentum 



k. With 

Ki = J l . 

we immediately find after performing the color trace 

wn(k, p) = 2 , w 22 (k, p) = 1 



wi 2 (k,p) = w 2 i(k,p) = 



(B5) 
(B6) 



2. Planar phase 



In the planar phase, we have 
1 



V + - 



Jf (1 - kf ) + Ji (1 - k z 2 ) - { Ji, J 2 } fcifc 2 + h fca7 7 7 ' k 



73+ — i _ 73+ 

' k,2 — 1 ' k,l 1 



(B7) 



where {Ji,J 2 } denotes the anticommutator of Ji and J 2 . In order to perform the color traces, we use Jf = Jf 
(i = 1,2) and 

Tr c [jf]=2 (i= 1,2,3), Tr c [J 2 jf] = 1 , Tr c [{ J u J 2 } 2 } = 2 , 
Tr c [J? J 3 ] = Tr c [j/{ J u J 2 }] =0 (* = 1, 2) . 

The only apparently additional Dirac trace which occurs in this phase is in fact identical to the above one, 

7 A (l-7 5 )7 5 7-kA k7 CT (l- 7 5 )7 5 7-pA p l = T ACT (k,p) . 



Tr 



D 



Consequently, with Eq. i|B4|) we find 



where 



?/(k, p) 



4fc 3 p 3 + (fciffi + fc 2 p 2 ) 2 - (fcip 2 - fc 2 pi) ; 



(B8) 
(B9) 



a; u (k,p) = i[3 + »j(k,p)] , w 22 (k,p) = i[l + 77(k,p)] , w 12 (k,p)=a; 2 i(k,p) = i[l-?7(k,p)] , (BIO) 



(Bll) 



3. j4 phase 



In this case, there arc three different quasiparticlc branches. Hence, we have three projection operators, 



1 



73+ 



J 2 , 



l + sgn(fc 3 )7 U 7 J 7-k 



vi 



1 



1 - sgn(fc 3 ) 7 U 7 5 7 ■ k 



k,2 2 - 

Again, all Dirac traces can be reduced to the previous one, because of Eq. (|B9|I and 



Ks = 1 - J: 



Tr 



7 A (l-7 b )7 5 7-kA k7 CT (l-7^ 7 u A r 
= Tr D [ 7 A (1 - 7 5 ) 7° A k 7 CT (1 - 7 5 ) 7 5 7 • P A p ] = T ACT (k, p) 
We thus find immediately, again using Eq. (|B4(1 . 



wn(k,p) = -[l + sgn(fc 3 )][l+sgn(p 3 )] , w 22 (k, p) = -[1 - sgn(fc 3 )][l - sgn(p 3 )] 
wi 2 (k,p) = -[l + sgn(fc 3 )][l -sgn(p 3 )] , w 2 i(k, p) = -[1 - sgn(fc 3 )][l + sgn(p 3 )] 
w 33 (k,p) = 1, Wi 3 (k, p) = W3i(k, p) = w 23 (k,p) = cj 32 (k,p) = . 



(B12) 



(B13) 

(B14) 

(B15) 
(B16) 



All vanishing functions Lu rs (k,p) are zero because of a vanishing color trace. This is exactly as in the polar phase 
where also all functions uj rs (k, p), that include one gapped and one ungapped branch, vanish. 
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4. CSL phase 



In the CSL phase, 



^ = 1+2 



J -7x0*0 



It is convenient to write these projectors with the help of color indices a, b = 1, 2, 3, 



_ = S ab + ^7-L,6(k)7_L,o(k) , 

30 Z 



L = -^7±,b(k)7±,a(k) 

3,0 Z 



(B17) 



(B18) 



Using these expressions, we compute the color and Dirac traces. The only additional nontrivial trace that occurs can 
again be expressed in terms of the previous one, 



J2 Tr c [7^1 -7 5 )7 7±,a(k)7±, 6 (k)A k7 ff (l-7 5 ) 7 7i , b (p)7 ± , Q (p)A p j = (1 + k • p) 2 T Xa (k, p) . (B19) 



.,6=1 



Then, with Eq. (|B4fl wc find 



wu(k,p) = i + i(i + k- P ) 2 , ^ 22 (k,p) = i(i + k- P ) 2 

w u (k,p) = w 21 (k,p) = l-i(l + k-p) 2 . 



Note that we have also used the following identity, 

3 



E K«(k) 



= -2 



(B20) 



(B21) 



APPENDIX C: ANGULAR INTEGRALS 



Here we compute the integral 

I(p„)= I dftp I dSl k (l - cos0„ d ) 5(cos8q - cosO ud ) F(k,p) 



(CI) 



The result shall be applied to Eq. Thus, the function F(k, p) stands for (k, p,u), defined in Eq. i.e., 

for notational convenience we omit all indices since only the dependence on k and p is of relevance. For the definition 
of cos 6*o see text below Eqs. (|31H . 

There is a fixed direction in position space, defined by the order parameter. We choose our frame such that the 
z-axis points into this direction. Moreover, for the dfl^ integration, we may choose the frame such that the three- 
momentum of the d quark, p, lies in the xz-pl&ne. Then, we can write the unit vectors of the momenta of u quark, d 
quark, and neutrino as 



sin t) u cos tp u 
sin 9 U sin ip u 
cos 6„ 




sin v v cos y) v 
sin 9 V sin ip v 
cos 6*„ 



(C2) 



respectively, and cos 6 u d = k ■ p, cos Q v d = Pv • P- (In this appendix, no confusion of the angle tp u with the ratio 
= 4>u/T is possible.) In the chosen frame, the most general angular dependence of the function F, which applies 
for all phases we consider, is 



F(k, p) = F(cos 9 U , cos 8d, sin 9 U sin 6d, cos tp u ) 



(C3) 



The angular dependencies enter this function through the eigenvalues Ak.r, A p . s and the functions w rs (k, p). In the 
CSL phase, the eigenvalues are constant. However, the functions uj rs depend on the scalar product k • p, wherefore 
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the arguments sui9 u sin6 d and cos <p u enter the function F. In all other phases, Ak, r and A PjS are angular dependent. 
They depend only on cos9 u and cos (9^ respectively. The functions uj rs are constant in the polar phase, whereas they 
depend on cos 9 U and cos 9 d in the A phase. In the planar phase, the angular dependence of the functions u> rs seems 
complicated, see Eqs. (|B10|) and (IB11|) . However, with f>2 = 0, the only additional argument is cosip u . 
Let us abbreviate Ci = cos^,; (i = u, d, 0). Then, 



^(Pi/) = J dn p (1 - cose vd ) J dip u J d^ u 5 (^o - C«Cd - ^(1 - £2)(1 _ Q cos ip u ^j 
Due to the ^-function, the integral is only nonzero if 

Co - CuCd 



va-axi-a) 



< i 



or, equivalently, if C u < £ u < £+, where 



Consequently, 



I(p„) = J / dSlp (l-coB0 vd ) J < 



t f y/(i - - e d ), Vi t~£\te d ) ) 



(C4) 



(C5) 



(C6) 



V(i-e2)(i-^)-(e -^ d ) 2 



(C7) 



In order to perform the C« integral explicitly, we introduce the new variable £' u = (£„ — £o£d)/y/0- — Co)(l — Cd) w ^ n 
the integration range from —1 to 1. In terms of the new variable, the integrand becomes a non-singular function 
of Co when Co - >■ 1- Because of the kinematics of the Urea processes, see text after Eqs. I|29|). the actual value of 
Co — 1 — K Me/(A t ii/ i d) is very close to 1. Therefore, to leading order, we may set Co = 1- Then, the integral over £' u 
can be performed analytically, and we arrive at 



/( Pl ,)~ 2tt J dVt p {l- C0S e vd )F{Uid,l-i 2 d ,l) 



(C8) 



The arguments in the function F show that the leading result is obtained from k ~ p. Applying this approximation 
to the results of the previous appendix, we see that the functions tu rs become constant in the case of the planar and 
CSL phase. (Note, in particular, that ?y(k,k) = 1 in the planar phase.) 

In order to perform the dQ p integration, we have to reinstall the second component of the d quark momentum, 
f>2 = sin#d siniyjd, in the term cos9 vd . We may perform the integral over (p d and obtain 



I( P „)~4tt 2 / dC(l-Ccos^)f(C,C,l-C 2 ,l) 



(C9) 



APPENDIX D: EMISSIVITY AT LOW TEMPERATURE 



In this appendix, we compute the tp dependence of the function G(<p,ip), cf. definition i|41|> . for ip — > oo. This 
corresponds to the behavior of the neutrino emissivity of the gapped modes for small temperatures. For the phases 
in which the gap function has no nodes (CSL and planar phases), it is useful to compute the asymptotic behavior of 
the following integral separately, 

POO POO POO , _} _1 

I(<p,<p) = V / dvv 3 / dx dy ( e - e ^y 2+l ? 2 + 1) (e e2 v * 2 +v 2 + 1) 

e u e 2 =± Jo Jo Jo ^ ' V ' 
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After performing the summation over e\ and e2 explicitly and taking into account that exp(y / a; 2 + p 2 ) 1 and 
exp(^y 2 + ip 2 ) 3> 1, we arrive at the following expression, 

poo />oc poo / 2 1 

I(ip,<p) ~ / (itin 3 / da; / e?y 



JO \ e V^+^ + e *> + vV+¥^ e -u+- v /a: 2 +v 2 _)_ g-y/l^+V 2 

1 1 \ 



3 v + e v5v+V»V g.+ViH^+V^+p 2 



(D2) 



Both the third and fourth terms in the integrand yield contributions of the order of exp(— 2tp). For the third term we 
obtain 



oo poo -i poo roc 



dvv* I dx I dy ■ — — — ■ = -6 / dx I dye-V* 2 +v 2 -Vy 2 +¥> 2 

Jo Jo e* + e V^+^+V^+^ Jo Jo 

x Li^-eV^+^+vV+v 2 ) 
~ 2irtp 5 e-' 2lp , (D3) 

where we used the asymptotic behavior of the poly logarithm Li4(— e z ) ~ — j^z 4 for z — > oo. The fourth term yields 

oo r oc i 



dvv / dx dy -j^=. — ~ 37ry exp(— 2(p) (D4) 

JO JO e v + ytx 2 +ip 2 + y/y 2 +ip 2 

We neglect both contributions (|D3|1 and (|D4|) to the integral I(ip, p) since they are suppressed stronger than the first 
two terms terms in Eq. I|D2(1 . The latter yield contributions of the order of exp(— <p), 

poo poo poo p~^P 

I(tp,<p) ~ 2 / dvv 3 dx dy 
Jo Jo Jo 

/>oo />OC fOO 

= *pe-*J o dvv 3 j Q dxj o dy ex2+ev+y2 
~ 42.55 tpe-f, (D5) 



e x 2 /(2 V ) + e v+y 2 /{2 V ) 
1 



where the overall number is obtained by performing the numerical integration. 



CSL phase 



Using Eqs. Q37|). I|41a|l . I|D5|I. and A^i = 2 (see Table |nj we immediately conclude 

5040 
457 7T 6 

0.69 tpe-^f . (D6) 



2. Planar phase 



For the planar phase, we may also use Eq. (|D5|I because the gap function has no nodes around the Fermi sphere. 
However, the result contains an additional angular integration. Using Eqs. Ij37|) . Q41a[) . and A^i = 1 + £ 2 , we derive 

2520 



G( ^ ) = 457^6 J ^IiV^ + e^V^+ev)- (D7) 



Approximating the angular integral as 



dt^T+ee-V^e* ~ / dee~ (1+ « 2/2)v = V^ — , (D8) 

J-oo 

we obtain 

G(ip,tp) ~ 0.61 ^pe-'f . (D9) 
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3. Polar phase 



In the polar phase, the low temperature behavior of the neutrino emissivity is dominated by the regions around the 
nodes of the gap function. Since A^i = sin 2 (9^, the dominant contribution comes from the vicinities of the north and 
south pole of the Fermi sphere. The gap behaves linear around these points, •v/Ak^i — To obtain the leading order 
result in tp, we integrate over the region where the quasiparticle energy is less than or of the order of the scale set 
by the temperature, i.e., ^/Ak,i </> < 7rT, or, equivalently, #k tt/V (and analogously for the south pole of the Fermi 
sphere). Consequently, we may use Eq. H41a|l . set tp u = ip,i = and integrate over the appropriate angular regions. 
Then, upon using ujh (£) = 2 (see Table HU and the integral J2B we obtain 

G(cp, tp) ~ / d8 k # k = — ~ — . (D10) 
Jo 2l P V 



4. A phase 



As in the polar phase, the dominant contribution comes from the nodes of the gap function. Therefore, we can 
neglect the term corresponding to the first (fully gapped) branch in Eq. I|41b|) . We keep the term that corresponds to 
the gap function given by Ak,2 = (1 — | cos#k|) 2 - In contrast to the polar phase, the gap function behaves quadratically 
in the angular directions in the vicinity of the nodes, y/\u.,2 — S 2 ./2. Therefore, in order to be consistent with the 
estimate in the polar phase, we restrict the integral to the regions < 8^ < y/2n/tp and it — y/2ir/ip < 8^ < tt. With 
^22(0 = 1 — s g n (C) we obtain 

1 r* 7t 1 ^7 

G((p, tp) ~ - d 8 k (n-8 k ) = — c— . (Dll) 

Note that due to the factor 0*22 m f ac t only one of the nodes, 8^ = tt, yields a nonzero contribution. The physical 
reason for this is explained in detail in Sec. IVII 



APPENDIX E: SPECIFIC HEAT AT LOW TEMPERATURE 



Here we compute the behavior of the function K(tp), cf. definition l|47|l . for tp — ► 00. Physically, this corresponds 
to the behavior of the specific heat of the gapped modes for small temperatures. From Eqs. I|47[) we conclude that for 
large tp and hence (p 2 3> </> 2 /T 2 

3 [°° f 1 1 

K(<p) - / dx / d ^ ^ + A «.! n 1 9 ; == ^\ ( CSL ' P lanar > P° lar ) > ( Ela ) 

2tt 2 Jo J-i 1 + cosh( v /z 2 + A^i^ 2 ) 

3 2 fOO p\ -1 

X ^ / 9 / , j. / „ 9\ J- 



*w - ggy„ ■ iai2+W) i + co 8h( ^ + A,,^ <Eib) 



1. CSL phase 



In the CSL phase, A^i = 2, thus the angular integration d£ is trivial. Moreover, we may approximate 

1 ~ 9 -y/x^+2^ ^ 9 -V2 



1 + COsh (y/x 2 + ^t,l<P 2 

Consequently, 



2 e rv x +2ip ~ 2e" v zlp ~^7^ . (E2) 



K(<p) ~ A f" dx ( x 2 + 2 2) e -V^^-^ „ 3v| ( ^ )5/2 e -V2o „ Lgl ^5/2 ^ _ (E3) 
7T J /r, 77" / 
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2. Planar phase 

As in the CSL phase, the gap has no nodes in the planar phase. Therefore, a similar approximation can be used, 

o />1 roc x 2 

K((f) ~ — / dS,e-V T +^'P / dx [x 2 + (1 +C 2 )^ 2 ] e~*Vi+(^ 
n " J-i Jo 



V2tt 3 /- 



5 / 2 / d{,e-^+?^{l+ef A , (E4) 



where the contribution originating from the term ~ x 2 in the square bracket of the first line of Eq. (|E4|) . has been 
neglected since it is suppressed by a power of (p. We can further approximate this expression to obtain 

K^c—^—ip 5 / 2 ^e- (1+ «" /2)v = -^e- v ~0.95(p 2 e-^ . (E5) 

V2 7T 3 / 2 J^oo IT 

3. Polar phase 



In the polar phase, with Ak.i = sin #k, we have 

= TTT / dx dd* sm9 k (x 2 + sin 2 # k v? 2 ) — . (E6) 

27r Jo Jo 1 + cosh yi 2 + sin 2 k <p 2 



In this case, the dominant contribution comes from the vicinities of the north and south pole of the Fermi sphere, 
where the gap vanishes. The gap behaves linear around these points, «/ Ak,i — #k- To obtain the leading order result 
in tp, we integrate over the region where the quasiparticle energy is less than or of the order of the scale set by the 
temperature, i.e., yAk^i^o wT, or, equivalently, #k J; ir/ip (north pole of the Fermi sphere, the south pole yields 
the same result). Consequently, 

4. phase 

In the A phase, we use Eq. i|Elb|) . i.e., in principle, we have to consider two different angular gap structures. 
However, the small temperature behavior is dominated by the nodes of the gap. Therefore, it is sufficient to consider 
only the second branch, Ak,2 = (1 — | cos6>k|) 2 , 



/>0O />7T 

/ dx d6 k sin0 k [x 2 + (1 - | cos0 k |)V ] - 
Jo Jo 1 



3 

4tt 2 Jo Jq j 1 + cosh (y/x 2 + (1 - | cos6» k |) 2 <p'' 



K{ip) = — j dx j d6» k sm(9 k [x 2 + (1 - | cos0k|)V ] — ■ ( E8 ) 



Contrary to the polar phase, the gap function behaves quadratically in the angular directions in the vicinity of the 
nodes, ^/Ak,2 — $k/2- Therefore, in order to be consistent with the estimate in the polar phase, we restrict the 
integral to the region # k ^ y/2ir/ip, 

Q roc 1 py/ln/p / ai \ _ 014 

K ^^l ^TT^ihwl *<>4 2 + M = ^V- (E9) 
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